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Abstract 

We evaluate the small-amplitude excitations of a spin-polarized vapour 
of Fermi atoms confined inside a harmonic trap. The dispersion law 
u — tof[l + 4n(n + I + 2)/3] 1//2 is obtained for the vapour in the colli- 
sional regime inside a spherical trap of frequency ujf, with n the number 
of radial nodes and I the orbital angular momentum. The low-energy ex- 
citations are also treated in the case of an axially symmetric harmonic 
confinement. The collisionless regime is discussed with main reference to 
a Landau-Boltzmann equation for the Wigner distribution function: this 
equation is solved within a variational approach allowing an account for 
non-linearities. A comparative discussion of the eigenmodes of oscillation 
for confined Fermi and Bose vapours is presented in an Appendix. 

Pacs:67.40 Db Quantum statistical theory; ground state, elementary excita- 
tions. 



1 Introduction 

The achievement of condensation in confined vapours of bosonic alkali atoms pj- 
[[| and more recently in vapours of Hydrogen |Q] has been giving great stimulus 
to the study of many-body and quantum statistical effects in dilute quantal flu- 
ids. In particular, the observation of small-amplitude shape-deformation modes 
in these bosonic condensates at very low temperature ^|-|| has provided a cru- 
cial test for the theory as developed by many authors (sec for istance 
Trapping of fermionic species has also been reported for 6 Li (To) and 40 K jyj. 
Cooling of these Fermi gases into the regime of quantal degeneracy still remains 
to be achieved experimentally, but one may already anticipate that the observa- 
tion of the dynamic response of the vapour of fermionic atoms to modulation of 
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the trap will provide an important method of diagnostics. It is to the theoretical 
study of the elementary excitations in such dilute normal Fermi fluids at very 
low temperature that the present work is addressed. 

We assume that the vapour has reached thermal equilibrium inside a mag- 
netic trap and regard it as a dilute, polarized Fermi fluid. The s-wave collision 
between pairs of atoms are inhibited by the Pauli princile, so that to leading 
order only p-wave scattering and dipole magnetic interactions are important. 
However, these interaction effects are very small at the temperatures of present 
interest jl2|]-[[l3). The confined vapour may then be treated as a non-interacting 
Fermi gas, with resistance to shape deformation being provided by exchange in 
the regime of quantal degeneracy. 

In Sect. 2 we start treating an inhomogeneous gas composed of (a large 
number of) non-interacting fcrmions from the general viewpoint of the equation 
of motion for the Wigner distribution function and of its projections yielding 
the conservation laws in the form of quantal hydrodynamic equations |Q. On 
neglecting temperature fluctuations, the kinetic stress tensor contains the essen- 
tial physical input for the study of both the equilibrium density profile and the 
dynamics of density fluctuations. We assume next that (i) the collision mech- 
anisms that have brought the vapour to thermal equilibrium are still operative 
during its dynamics, and (ii) the inhomogeneity is weak on the lenght scale 
set by the inverse Fermi wavenumber. Under these approximations the kinetic 
stress tensor is simply related to the local Pauli pressure and to the velocity 
field. We recover from this approach the Thomas-Fermi form of the equilibrium 
density profile already found by M0lmer |l5| ] for fermions in a boson-fermion 
mixture at zero temperature. We proceed from here in Sect. 3 to evaluate the 
eigenmodes of small-amplitude oscillations for the Fermi vapour at zero temper- 
ature in a spherically symmetric trap. In an Appendix we show how one may 
recover within the same approach the results of Stringari |jj for the collective 
excitations of a Bose condensate in the strong coupling regime and compara- 
tively discuss the nature of the eigenmodes solutions for the Fermi gas and for 
the Bose condensate. 

In Sect. 4 we extend the discussion of the small-amplitude vibrations of an 
inhomogeneous Fermi gas to the case of confinement within an axially sym- 
metric harmonic trap and give explicit results for the low-lying eigenmodes. In 
Sect. 5 we consider instead the dynamics of a Fermi gas under spherical con- 
finement in the case where collisions have a negligible effect on the deviations 
from the equilibrium. We make brief reference to zero-sound excitation of the 
homogeneous Fermi gas in such a collisionless regime and then focus on a vari- 
ational treatment of a Landau-Boltzmann equation for the Wigner distribution 
function. Finally, Sect. 6 gives a brief summary of our main results and offers 
some concluding remarks. 
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2 Small amplitude excitations in the collisional 
regime 

We describe the inhomogeneous spin-polarized fermion cloud by the one-particle 
density matrix 

p(x,x / ;t) = (*t( X)t )*( x ',t)) (1) 

where if? and if!' are the fermionic field operators satisfying the usual anti- 
commutation relations, and we have suppressed the spin index, which is frozen 
because the magnetically-trapped fermions belong to a single hyperfine level. 
The equation of motion for the density matrix of non-interacting fermions in an 
external potential V ex t(x, t) is (see for istance Singwi and Tosi Jig]): 

-2 -| 

p(x,x';t)=0 (2) 



ihd t - — (vi - vl,) + v ext {-K, t) - v ext { x ', t) 



We shall take V ex t(x, t) as the sum of a static harmonic potential Vb(x) and the 
perturbing potential V% (x, t) . 

The generalized hydrodynamic equations expressing the conservation laws 
are obtained by expanding Eqn. (||) in powers of r = x — x' around the diagonal 
R = (x + x')/2. We introduce the Wigner distribution function which is related 
to p(x, x', t) by: 

/ P (R, t) = J drp(R + r/2, R - r/2; t) exp(ip • r/H), (3) 
and its successive moments giving the density of particles n(R, <), 

n(n,t) = J ^3/ p (R,f), (4) 
the current density j(R, t), 

and the kinetic stress tensor Ry (R, t) , 

Then by the indicated expansion procedure we obtain the conservation laws in 
the form: 

d t n(R,t) = -V ■ j(K,t) (7) 

and 

md t j i (R,t), = -V i n i . / (R,t) -n(R,t)ViF ext (R,t) . (8) 

We consider below only isothermal fluctuations and therefore omit to write 
down the energy conservation equation. 
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2.1 Approximations on the kinetic stress tensor 

The kinetic stress tensor of an homogeneous fluid is related to the pressure P 
and to the velocity field v by 

Hij = PSij + nmViVj (9) 

Fore the ideal spin-polarized Fermi gas at zero temperature we have 

P = 2E/W = 3ne//5 (10) 

where E is the gound state energy and e/ = h 2 (6ir 2 n) 2 / 3 /2m is the Fermi energy. 
On the other hand, the second term on the RHS of Eqn. (||) can be dropped in 
the linear regime. 

These results can be used to obtain the kinetic stress tensor of the inhomo- 
geneous Fermi fluid if the inhomogeneity is sufficiently weak to allow a local- 
density approximation, i.e. 

TLij{R,t) = 8 l3 P{R,t) = S^A[n(R,t)} 5 / 3 (11) 

where A — fi 2 (dm 2 ) 2 / 3 /2m. More explicitly, we are assuming that the lenght 
scale for the variation of the density profile in space is large relative to the 
inverse Fermi wavenumber kj ■ 

We remark that insertion of Eqn. ( |ll| ) into Eqns. (0) and (||) and reduction 
to the case of a static external potential yields the equation for the equilibrium 
density profile rio(R) in the form 

-AVnJ 3 + n Wo = 0. (12) 
o 

The solution of Eqn. ([l2|) is the Thomas-Fermi profile for the confined Fermi 
vapour: 

3/2 



f-Vo(R) 



n (R) = 9(e f V (R)) ^ 1 ^ ' j (13) 

(see for instance M0lmer @). The profile vanishes continuosly at a radius Rf, 
say. 

2.2 Equation of motion for small amplitude fluctuations 

With regard to the dynamics of the fermionic cloud, the approximation ( p"2| ) 
is most suitable for treating a "collisional" regime in which the relaxation of 
the fluctuations is rapid on the time scale set by the frequency of the driving 
potential. We shall have to return to this point in Sect. 5 below. Here we 
proceed to study the dynamics of the fermion cloud for small distorsions of its 
density profile around the equilibrium Thomas-Fermi profile given in Eqn. ( |l3|) . 

If we denote by ni(R, t) such a small-amplitude density fluctuation, Eqns. (|7j)- 
(||) and ( pi] ) can be linearized and combined to yield the equation of motion 

dfm(R,i) = -V R • (ni(R)V R y (R) + ^V R L 2/3 (R)m(R)l } . (14) 
ml 3 L J I 
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Equation (Q) holds at resonance, since we have dropped the external pertur- 
bation drive. In the homogeneous limit (Vb(R) — > 0) it describes the sound 
propagation at a speed c 2 = 2e//3m - a well-known result for first sound in the 
ideal Fermi gas fl?]]. In the so-called collisional regime the local compression 
and expansion of the homogeneous fluid which are induced by the sound wave 
arc reflected in spherically symmetric oscillations of the diameter of the Fermi 
sphere. 



3 Small oscillations in a spherical trap 



We determine in this section the solution of Eqn. fll4| ) for the case of a spher- 
ically symmetric harmonic confinement, i.e. Vo(R) = muj 2 R 2 /2. We Fourier- 
transform Eqn. (|lj) with respect to time and rescale the lengths with the har- 
monic oscillator length a^ = (ft/mwy) 1 / 2 , by setting R = xa/jo. We get: 

- 3w 2 ni(x,w) = cu 2 f \7 x • [ni(x,w)V x (> 2 /2) + (X 2 - .x 2 )V x ni(x, to)] (15) 

where X — ^/2ei?/mw 2 /a/ lo is the scaled radius of the fermionic cloud. We 
search for spherical solutions of Eqn. (|l^) having the form: 

m(x, u) = 9(X - x)x l F(x/X)Y l m (9, 4>) (16) 

with Y, m (9, tf>) being the spherical harmonics. Clearly, we have imposed on the 
eigenmodes solutions the condition that they be non-vanishing only inside the 
equilibrium cloud radius X. In addition, we shall impose that the solutions 
vanish continuously at the cloud boundary x — X, as a consequence of Fermi 
statistics giving a high cost in kinetic energy to rapid variations of the density 
profile in space. This boundary condition on the eigenmodes of vibration for 
a trapped Fermi gas is further discussed in the Appendix and it is contrasted 
there with that imposed e. g. by Stringari on those of a trapped (interacting) 
Bose gas, where considerations of kinetic energy only play a secondary role. 

After setting y = x/X and substituting expression (|l6|) into Eqn. (|l5| ) we 
find the differential equation satisfied by F(y): 

h ^ d 2 f 2(l + l)-(2l + 3)y 2 \ d 

+ (3 + H"/u Jf ) 2 -l)F(y) = 0. (17) 

with the boundary condition F (1) = 0. We adopt the Fuchs method Q for 
solving in a series form ordinary differential equations around regular singular 
points. Using the above-mentioned condition of continuity at the boundary, and 
the fact that F{y) is an even function of its argument, the solutions of Eqn. ( p"7| ) 
are 

oo 

F(y) = (l-y 2 ) 1 ? 2 Y / Ml-y 2 ) k (18) 

k=0 
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and the coefficients au satisfy the recurrence equation: 
(2k + l)(2fc + 21 + 3) - (3 + 3(oj/uj f ) 2 - I) 
ak+1 = 2(2fc + 3)(A: + l) ak (19) 

Finally, the dispersion relation for the normal modes is obtained by imposing 
that the solutions (|l8|)-(|l9|) reduce to polynomials, i.e. a n +i = for an integer 
n representing the number of internal nodes of the density fluctuation profile. 
This yields 

f 4 \ 1/2 
u(n,l)=Uf U + -n{n + l + 2)\ (20) 



The dispersion relation (|2C|) correctly displays the sloshing mode solution (n = 0, 
I = 1, m = 0) at a frequency to = u)f, in agreement with the generalized Kohn 
theorem pj| . 

It is also interesting to notice that the eigenfrequencies surface modes at 
n = are the same as those found by Stringari JtJ for a Bose condensate in 
a spherical trap and that the I — modes resemble those found for a Fermi 
superfluid by Baranov and Petrov [M . 



4 Extension to axially-symmetric magnetic con- 
finement 

We evaluate in this section the low-lying solutions of Eqn. (|l4|) in the case 
of axially symmetric harmonic confinement, which is more directly relevant in 
regard to experiment. The confining potential is chosen as 

V ( Pl z) = ^mcu 2 f (p 2 + X 2 z 2 ) (21) 

where p = \J x 2 + y 2 is the radial coordinate and A is the anisotropy parameter. 

In the axially symmetric case only the L z component of the angular momen- 
tum along the z direction is conserved, whereas I is no longer a good quantum 
number. However, for the low-lying modes we can still start from a spherical 
base in Eqn. (|l^), and suitably modify it so as to account for the cylindrical 
symmetry of the problem. As we shall see through an example below, the con- 
sequence of the reduced symmetry is a coupling between some of the modes 
which are characterized by different values of I in the spherical base. 

We give here the explicit results for the eigenfrequencies and the density- 
fluctuation profiles in some istances of dipolar and quadrupolar modes In the ex- 
pressions reported below, all lengths are still scaled in units of ato = [h/mujf) 1 ^ 2 . 

4.1 Dipole sloshing mode in the z direction. 

Starting from the spherical-base solution corresponding to (n = 0, 1 = 1, m = 0), 
we construct the density profile 

n^\p,z) = {l- P 2 -X 2 z 2 f/ 2 P (22) 
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The corresponding frequency is io 2 = A 2 w 2 in agreement with the generalized 
Kohn theorem jl9| . 



4.2 Quadrupole surface modes. 

The anisotropy introduces a splitting between the quadrupolar modes corre- 
sponding to (n = 0, 1 = 2, m = 1) and to (n = 0, 1 = 2, to = 2). the density 
profiles and the eigenfrequencies are as follows: 



if' 2 ' 1} (p,z) = (l-p 2 - \ 2 z 2 ) x / 2 pzvq>{i<i>) (23) 
lu = (l + A 2 ) 1/2 w/ (24) 



r4°> 2 < 2 )(p, z ) = (1 - p 2 - \ 2 z 2 fl 2 p 2 exp(2*0) (25) 

uj = 2 1/2 uj f (26) 

In these equations <f> is the angular momentum variable in the transverse plain. 
The two eigenfrequencies become degenerate in the spherical limit. 

4.3 Monopole and quadrupole coupled modes. 

The anisotropy introduces a coupling between the modes described by the quan- 
tum numbers (n = 0, 1 = 2,m = 0) and [n = 1, 1 = 0, m = 0). The associated 
density fluctuations are 

„(°A°) = (1 - p 2 - A 2 z 2 )!/ 2 (2z 2 - p 2 ) (27) 

and 



(1,0,0) 



(l-p 2 -A 2 z 2 ) 1 / 2 (a-p 2 -z 2 ) (28) 



where a has to be determined by a normalization condition. For the coupled 
eigenmode equations we find the eigenfrequencies as: 



2 2 
UJ = LO 



4A 2 + 5 ± V16A 4 - 32A 2 + 25 



(29) 



The above examples should be sufficiently illustrative of the procedure to be 
followed in determining the small-amplitude oscillations in the case of axially 
anisotropy. 
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5 Collisionless regime 



We return to the case of spherical confinement for the purpose of discussing 
the role of collisions in the dynamics of the fermion cloud. We have already 
remarked in Sect. 2 that in the homogeneous limit Eqn. @, which was derived 
with the help of Eqn. ( |ll|) for the kinetic stress tensor, describes first sound 
in a collisional regime where the local Fermi sphere executes a "breathing" 
oscillation. In the so-called collisionless (or zero-sound) regime, on the other 
hand, the relaxation time of density fluctuations in the homogeneous Fermi 
fluid is long compared with the period of the eigenmode and the local Fermi 
sphere undergoes anisotropic oscillations (see for istance Pines and Nozieres 
[p7r| ) . In the limit of vanishing coupling the speed of the zero sound tends to the 
Fermi velocity vf (compared with the value 3~ 1 ^ 2 vf for the velocity of the first 
sound) and the deformation of the Fermi surface reduces to a small bump in the 
direction of the propagation vector. The collective zero-sound mode involves in 
this limit only a small number of quasiparticles and essentially propagates at 
their velocity. 

For a discussion of the dynamics of a magnetically confined Fermi cloud in 
a regime where collisions have little influence, we resort to an equation for the 
Wigner distribution function having the form of a Boltzmann equation with the 
collision integral set to zero, that is 

d t f p (R, + - ■ V R / P (R, f) - V R Vo(R) • V P / P (R, t) = (30) 

TO 

This equation is in fact analogous to the Landau transport equation for quasi- 
particles in a homogeneous Fermi fluid in the collisionless regime and is expected 
to be valid down to very low temperature provided that the cloud contains a 
large number of fermions. 

In searching for solutions of equation ( |30|) , we focus on the shape-deformation 
modes which can be monitored by measuring the mean square radius of the 
cloud as a function of time. Within this class of modes, we adopt a variational 
Ansatz as already proposed by Bijlsma and Stoof [^3], |24| for a Bose cloud. This 
amounts to choosing / P (R, t) in the form of the equilibrium solution $ eg (p, R) 
after rescaling both the space and momentum variable through time-dependent 
factors in each geometrical direction. More precisely, for (i — x, y, z) we intro- 
duce variational parameters at(t) describing the deviation of the mean square 
sizes of the cloud from their equilibrium values, 

Oi{t) = [{Ik{tY)/{R*) eq ]V* (31) 

The form of / P (R, t) is taken as 

/ P (R, t) = (c x CyC z )$ eq (sJciotiipi - mRidi/ai); ^/c^Ri/ai) . (32) 

the time-independent coefficients c; being introduced in Eqn. ( |32| ) in order to 
ensure that the initial momentum distribution is isotropic. The equilibrium 
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form of the Wigner function for the degenerate Fermi gas subject to spherical 
confinement can in turn be chosen as a generalized Fermi sphere: 

<^ eq (p,R)=e(j^ + ^mu; 2 R 2 -^ (33) 



/if being the chemical potential. Of course, Eqn. ( p3|) yields the same equilib- 
rium densitty as that reported earlier in eq. (13|). We emphasize that in the 
present treatment we are assuming that, even though collisions have negligible 
influence on the dynamics of the cloud, they must have been active in the past 
in order to bring the cloud to the state of thermal equilibrium. 

The equations of motion for the aj(t) in each spatial direction are obtained 
by taking the moments of equation (|30|)witli respect to RiPi ||^, [2~|| : 

Ul 2 f 

di +oj}ai = (34) 
It is worth noticing that for a non-interacting gas Eqn. (p4) follows directly 



from the dynamical scaling ansatz (32) and from the virial theorem, and does 



not require a specific assumption on the shape of the variational function. 

The solutions of the non-linear equation of motion ([m]), under the initial 
condition ai(0) = 1, is 

Oi(t) = -J—[c + (c 2 - Aujj) 1 / 2 sm(±2to ft - 0)] 1 / 2 (35) 

with sin(</>) = [(c — 2uj 2 )/{c + 2a; 2 )] 1 / 2 and the constant c determined by the 
initial velocity according to: 

d(0) = ±(c-2lj}) 1 / 2 (36) 

A more direct contact can be made with the results reported in the preceding 
sections for the gas in the collisional regime by linearizing Eqn. (p4]). In such 
small-oscillation regime, Eqn. (Q) becomes an equation for the linearized mean 
square radius a^\t): 



n 



(1) I A..»J?) 



4Lojay> = . (37) 



Evidently one obtains an oscillation of the cloud at frequency w = 2a;/, corre- 
sponding to the frequency of a single-particle excitation in the cloud. 

As a final remark we point out that while Eqn. (^o|) satisfies the general- 
ized Kohn theorem, a more sophisticated variational ansatz would be needed to 
extract from it the sloshing mode urn . 



9 



6 Summary and concluding remarks 



In summary, we have presented an investigation of the collective excitations 
of a spin-polarized Fermi gas confined in an external harmonic trap. We have 
obtained analytic results for the small-amplitude eigenmodes in a dynamical 
regime corresponding to first sound in the bulk and proposed a Landau-Boltzmann 
approach to treat the dynamics of the confined gas in a dynamical regime cor- 
responding to zero sound in the bulk. In spite of the assumed negligible role 
of interactions between the particles of the gas in contributing to a mean field 
sustaining a collective oscillation, we have seen that the Pauli pressure alone 
can determine a rather rich dynamical behaviour. 

Trapped mixtures of bosonic and fermionic species are expected to become 
accessible to experiment in near furure. The investigation of collective excita- 
tions in such mixtures seems a natural direction of development for the present 
study. The density profiles of the mixtures have been determined theoretically 
both in the ground state at zero temperature [T^ | and in the equilibrium state 
at finite temperature JllJ . From these studies one can easily envisage situations 
in which the dynamics of a boson-fermion mixture should be amenable to ana- 
lytical approaches. Theoretical progress in this area will be reported elsewhere. 
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Appendix A. Collective excitations of Fermi ver- 
sus Bose-condensed clouds 



In this Appendix we first show how the dynamics of a dilute Bose condensed 
cloud at zero temperature can be studied within the Wigner-distribution for- 
malism that we have used for the spin-polarized Fermi gas (see also ju| ) . Wc 
thereby derive the analogue of the equation of motion (|14j ) for density fluctua- 
tions in the Bose-condensed cloud and point out that for harmonic confinement 
it coincides with the one obtained by Stringari in a hydrodynamic approach 
to the strong coupling limit. We then comparatively discuss the nature of the 
solutions of the two equations of motion for the Fermi gas and for the Bose- 
condensed cloud. 



Al. Equation of motion for density fluctuations in a Bose- 
condensed cloud 

In the case of a Bose-condensed system at zero temperature, and neglecting 
depletion due to the interactions, the one-body density matrix introduced in 
Eqn. (|l|) reduces to the simple product of two condensate wave functions, 

(*t(x, i)*(x', t)) = **(x, t)$(x' ) t). (38) 

The generalized hydrodynamic equations (0) and (|^) are thereby closed and 
the kinetic stress tensor is explicitly given by 



IXy(R,t) 



ft 2 



d d v v 

__*-(R-- >t )*(R + -, t ) 



(39) 

r=0 



However, for a Bose condensate with repulsive interactions in the strong 
coupling limit the kinetic stress tensor becomes in fact negligible compared 
with the confinement and interaction energy terms. The linearized equation of 
motion for density fluctuations becomes p5l 

fl£m(R,i) = -^J dx' J dx'> (R)V R( 5(R-x') 

+V R n (R)- VR<KR-x')Mx',x") m (x",i) ( 40 ) 
where i/(x',x") is the interparticle potential. 

Assuming contact interactions, i.e. f(x', x") = <7<5(x'— x") with g — 4irh 2 a/m 
where a is the s-wave scattering length, Eqn. ( (40|) immediately yields 

d t V(R,i) = -2-Vr • [no(R)V R ni(R,t)] (41) 
m 

This yields the equation of motion derived by Stringari [Q when the Thomas- 
Fermi solution is used for the equilibrium density profile hq (R) . Before proceed- 
ing we wish to emphasize that the nature of the Thomas-Fermi approximation 
is very different in the two cases that we are considering. In the spin-polarized 
Fermi gas the interactions are negligible, so that the kinetic term provides the 
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only energy scale. On the other hand, in a mesoscopic Bose condensate the 
interactions are dominant even at zero temperature. Nevertheless, the form of 
Eqn. (fH]) is essentially similar to that of Eqn. (fuj). 



A2. Comparative discussion of the solutions of Eqn. ( ^T| ) 
and Eqn. 

The solutions of Eqn. (^d]) for the Bose-condensed cloud vanish outside the cloud 
radius Rb and present a discontinuity at Rb . Indeed, by the Fuchs method 
used in Section 3 one may prove that it is not possible to impose continuity 
at Rb- The discontinuity is physically acceptable in view of the fact that the 
kinetic energy term has been set as negligible in taking the strong-coupling 
limit. On the other hand, Eqn. ( |l4| ) for the Fermion cloud admits, in addition 
to the solution that we have presented in the main text, a non- vanishing solution 
outside the cloud radius Rf. This solution is 

n(R,u) = C8(y - l)y-P+W u t)') . ( 42 ) 

Evidently, such a finite external solution going to zero only for y i — > oo as 



in Eqn. (42) is inconsistent with the assumption of small-amplitude oscillations 
of the cloud. Nevertheless, if for a moment one admits it and asks that the 
constant C be finite and determined by imposing continuity at the boundary 
with the internal solution (as necessary in a system where the kinetic energy is 
dominant), one finds a dispersion relation which is different from that given in 



Eqn. (gCj), that is 

3(1 +w 2 /^/) = l + An{n + l + l). (43) 

It is evident that the modes with n — and I < 3 are suppressed, and 
in particular the sloshing mode which is expected from the generalized Kohn 
theorem is absent. We have for this reason discarded this second solution and, by 
imposing that the internal solution vanishes continuously at the boundary (i.e. 
C = ), found the dispersion relation reported in Eqn. ( |2(i| ) in the main text. 
We also remark that each physically acceptable eigenfrequency of the cloud, 
with increasing n at fixed I in Eqn. (pC)), lies between two of the unphysical 
frequencies given by Eqn. ((43|). In conclusion, the different roles played by the 
kinetic energy in the Bose condensate and in the Fermi gas lead to very different 
forms of the density fluctuations in these two systems. In the former system the 
amplitude of the fluctuations vanishes discontinuously at the boundary of the 
cloud in the strong-coupling limit, whereas in the latter it vanishes continuously 
at the boundary. 
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